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Abstract 

The continuous infall of dark matter with low velocity dispersion from all directions in a galactic 
halo leads to the formation of caustics which are very small scale (~parsec) high density structures. 
If the dark matter is made up of SUSY neutralinos, the annihilation of these particles produces a 
characteristic spectrum of gamma rays which in principle could be detected. The annihilation signal 
at different energy bands is computed and compared with the expected gamma ray background. 
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I. INTRODUCTION 



According to the standard theory of cosmological structure formation, most of the mass in 
galactic halos is composed of non-baryonic matter with negligible primordial velocity disper- 
sion, commonly called cold dark matter (CDM), the constitution of which remains a mystery. 
Many theories have been proposed to account for the dark matter. Dark matter candidates 
may be non-thermal relics or thermal relics, depending upon the underlying particle physics. 
Non-thermal relics are particles which decoupled from the standard model particles without 
being in thermal equilibrium with them, the best example being the axion [lj]. Thermal 
relics on the other hand, were once in thermal equilibrium with other particles, but decou- 
pled when their interaction rate became comparable to the Hubble rate. These particles are 
collectively known as Weakly Interacting Massive Particles or WIMPs 0. 

There are many collaborations currently trying to detect dark matter, which include 
AD MX [3], DAMA/Nal [4], DAMA/LIBRA @, CDMS @, XENON Q, EDELWEISS 0, 
ZEPLIN [9|, etc. Dark matter detection experiments are usually of two kinds - direct 
detection experiments look for the recoil that occurs when a WIMP scatters off a target 
nucleus, while indirect detection experiments look for standard model particles which result 
from dark matter particle annihilation. In this article, we investigate the flux of gamma ray 
photons produced by WIMP annihilation in dark matter caustics. Previous work on particle 
annihilation in caustics includes 10, 11. 12l. 
by gravitational lensing 
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For the possibility of detecting caustics 
Caustics and their associated cold flows are also 
relevant to direct detection experiments [18, 19,(23, 21, 22 ■ 

Cold dark matter particles exist on a thin 3-dimensional hypersurface in phase space. To 
obtain the density in physical space, we must make a mappin g fro m phase s p ace to phys ical 
space. Caustics are locations where this mapping is singular [23L 24, [25|, 26, [27|, Ha 29, [3(J 

Caustics are made up of sections of the elementary catastrophes, 
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which are structurally stable and have a definite geometry that depends on the structure of 
the phase space manifold. In the limit of zero velocity dispersion, caustics are singularities 



36 



in physical space. In real galactic halos, the velocity dispersion cuts off the divergence 

It is important to ask whether the finite velocity dispersion of dark matter can soften the 
density contrast of caustics to the extent that the caustics are made irrelevant. The inner 
regions of phase space may well have thermalized over the age of the Universe. The infall 
of thermal dark matter does not produce observable caustics. However, it is believed that 
in large galaxies like our own, there is a lot of dark matter well outside the virial radius, 
yet gravitationally bound to the halo. These dark matter particles are in the process of 
thermalizing and have velocit y d ispersions that are small compared to the virial velocity 
dispersion of the galaxy j3l|, 33|. The dark matter particles reaching us today for the 
first time for example, are cold and are likely to produce a physically significant caustic. 
Similarly the particles that have fallen into the inner regions of the halo only a few times 
in the past may be expected to produce observable caustics. Thus caustics are formed if 
the outer regions of phase space are well resolved. An analysis of the rotation curves of 
32 galaxies 37j seems to provide evidence for the existence of dark matter caustics. The 
rotation curve of our galaxy also shows rises at the expected locations of the caustics 38 . 



Perhaps the best example of caustics on galactic scales is the occurrence of shells around 
giant elliptical galaxies [H, 0, 4l[ (e.g. NGC 3923, see |4l|). These shells are caustics in 
the distribution of starlight. They form when a dwarf galaxy falls into the gravitational 
potential of a much larger galaxy and is assimilated by it. The velocity dispersions of the 
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FIG. 1: Dark matter trajectories forming a caustic. The caustic is the envelope of the family of 
dark matter trajectories. Thus the dark matter density is large at the location of the caustic. 

stars of the dwarf galaxy are much smaller than the virial velocity dispersion of the giant 
galaxy. The infall of stars is therefore cold and presumably, collisionless. The occurrence 
of these caustics show us that cold flows are not unstable and that caustics do not require 
special initial conditions to form [33[ . This leads us to believe that the infall of non-thermal 
dark matter also produces caustics. 

There are two kinds of caustics - outer and inner 1321. The outer caustics are thin 



topological spheres surrounding galaxies (like the caustic of stars discussed previously). 
They are simple fold catastrophes and typically occur on scales of 100's of kpc for a galaxy 
like our own. The inner caustics have a more complicated geometry and are made up of 
sections of the higher order catastrophes [13]. Inner caustics typically occur on scales of 10 's 
of kpc. 

To see the formation of inner caustics, let us consider the infall of a perfectly cold flow 
of dark matter. If the infall is exactly spherically symmetric, the particle trajectories are 
radial and the infall produces a singularity at the center. If instead, the dark matter particles 
possess some distribution of angular momentum with respect to the halo center, the particle 
trajectories are non-radial, particularly in the inner regions of the halo. Fig. [1] shows an 
example of non-radial infall (in cross section) with the result that the dark matter density 
is enhanced along two thin fold catastrophe lines (surfaces in 3-dim. space) which meet at 



a cusp catastrophe point (line in 3-dim. space) [42|, |43|]. The caustic divides the xy space 
into two regions - one region with 1 particle trajectory passing through each point (outer 
region) and the other with 3 particle trajectories passing through each point (inner region). 
Close to the cusp, the dark matter trajectories are specified by the family of curves 

Vo x 

where (x = 0, y = 0) are the co-ordinates of the cusp and x , y are constants, a parametrizes 
the family. The caustic is the envelope of the family of curves, i.e. the locus of points tangent 
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to the curves. The envelope is obtained by solving Eq. [T] together with the tangency condition 
3« 2 



— = 0, using which we obtain the caustic curve 
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for x < 
0, we have 



The density is proportional to the sum J2i |3oj 2 ~ ( x / x o)\ where Oj are the real roots of 
the cubic Eq. [TJ When y — 0, the dark matter density is given by p(x, 0) = p c 

and p(x, 0) = p c + + §i) = 2 p c y for x > where p c is a constant. For x 

p(0, y) = — ^ . Close to one of the fold lines and in the inner region, the density falls 
off as the inverse square root of the distance to the fold along the direction perpendicular 
to the fold. For arbitrary (x,y), we need to solve Eq. [I] to obtain a%. The general solution 
is worked out in [loT ]. 

For arbitrary initial conditions, other catastrophes can occur. In general, we may 
parametrize the flow of cold dark matter particles by a three parameter label q(qj , a.2, 0.3) 
where 01,02,03 are three conveniently chosen parameters that describe the flow [32J. Let 
x(t, a) be the physical space location at time t, of the particle labeled a and is obtained by 
solving the equations of motion for given initial conditions. The physical space density is 
proportional to the sum 5Z i \dxfddi\~ 1 ( the Jacobian determinant of the map from a space 
to x space ) and is infinite (for a perfectly cold flow) at those points where dx/ddi = 0. For 
a perfectly cold flow, it can be shown that such a singularity necessarily occurs [34 ] . 

A perfectly cold flow cannot be realized in practice. The primordial velocity dispersion 
of WIMPs is estimated to be 5v = 3 x 10 _7 ^100 GeV/m x km/s Q. This causes the 
caustic surface to spread by an amount 5a ~ Sv x a/v (H| where a is the outer turnaround 
radius (few hundred kpc) of the flow of particles forming the caustic and v is the speed of 
the particles at the location of the (inner) caustic (few x 100 km/s (2l[). This minimum 
spread in the caustic surface sets an upper limit to the density. For the nearby caustics, we 
can estimate the spread in the caustic location 
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Let us now examine the effect of averaging the density over a finite region of space. 



Fig 2(a) shows the density along the line of sight, close to one of the fold lines of Fig [lj 



Close to the fold, but outside it, the dark matter density is very small. The density changes 
abruptly at the location of the fold and becomes very large (100 GeV/cm 3 was chosen 
as the density cut-off), decreasing smoothly thereafter. Figures 2(b), 2(c) and 2(d) show 



the density averaged over a cube of side 0.01 parsec, 0.1 parsec and 1 parsec respectively. 



In Fig. 2(b) , we see that the averaged density follows the true density faithfully, though 



the maximum averaged density does not quite reach the cut-off value. Fig. 2(c) no longer 
shows a sharp rise in density and Fig. 2(d) misses the caustic completely. We conclude that 
caustics are sub-parsec scale structures and are therefore difficult to resolve with large scale 
cosmological simulations which typically have spatial resolutions of 100's of pc. 
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FIG. 2: (color online) The effect of averaging the density over a finite volume. Fig. 2(a) shows the 
variation of dark matter density (in Gev/cm 3 ) with distance to the caustic ( in parsec ). A cut-off 
density of 100 Gev/cm was assumed. Figures |2(b)| |2(c)| and |2(d)| show the effect of averaging the 
density over a cube of side 5 for 5 = 0.01,0.1 and 1 parsec respectively. 



II. THE ANNIHILATION FLUX 



In the minimal supersymmetric extension of the standard model (MSSM), a good can- 
didate for the WIMP is the lightest neutralino which is a linear combination of the super- 
symmetric partners of the neutral electroweak gauge bosons and the neutral Higgs bosons. 
The characteristics of the annihilation signal depend both on the composition of the WIMP 
and its mass m x . The line emission signal ( XX ~~ > 77>XX ~^ %1 ) ls 1°°P suppressed and 
is therefore smaller than the continuum signal. The continuum flux (number of photons 
received with energies ranging from E\ to E 2 per unit detector area, per unit solid angle, 
per unit time) is given by [12|, |45|, |46[ 

$(E 1 ,E 2 ,6,<l>,An) = S{E u E 2 )x < L1 ' 
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dN^^/dE is the number of photons produced per annihilation channel per unit energy, b{ 
is the branching fraction of channel % and < av > is the thermally averaged cross section 
times the relative velocity. The factor of 2 in the denominator accounts for the fact that two 
WIMPs disappear per annihilation. The quantity EM(8, <fi) is called the emission measure 
and is the dark matter density squared, integrated along the line of sight, i.e. 

EM(6,<f>) = f dxp 2 (x). (6) 

J los 

and < EM > (6, 0, AQ) represents the emission measure in the direction (8, 0) averaged 
over a cone of angular extent AQ. We note that S depends solely on the particle physics 
while < EM > depends solely on the dark matter distribution. 



A. Estimating S 

Assuming that all the dark matter is composed of neutralinos, the quantity < av > is 
constrained by the known dark matter abundance [2j 



3 x 1(T 27 cm 3 s - 1 
< av > ~ — 



n x h 2 

-26 „„,3 „ -1 



3 x 10~ ib cm 3 s - 1 (7) 



The quantity dN^(Ej) / dx for the dominant channels may be approximated by the form [4^ 
48 . 4§] dN 1 /dx = ae~ bx /x where x is the dimensionless quantity E 1 /m x and (a, b) are 
constants for a g iven annihilation channel. The values of (a, b) for the important channels 
are given in [48(. Using these values, we may calculate the number of photons produced 
per annihilation within a specified energy range. Let us consider four energy bands: Energy 
Band I with photon energies from 30 MeV to 100 MeV, Band II with energies from 100 MeV 
to 1 GeV, Band III with energies from 1 GeV to 10 GeV and Band IV containing photon 
energies 10 GeV upto m x . The values of N 7 /m x are tabulated for the different energy bands, 
for m x = 50, 100, 200 GeV. 

TABLE I: N^/m^ in units of 10~ 4 GeV -2 for m x = 50 GeV 



Channel I II III IV 

WW, ZZ(0.73, 7.76) 106.8 85.2 18 0.52 
65(1,10.7) 146 114.4 21.6 0.32 

1,15.1) 160 122.4 19.6 0.12 

mt(0.95,6.5) 139.2 111.6 25.2 0.96 



B. Estimating < EM > 

The geometry of dark matter caustics depends on the spatial dark matter velocity distri- 
bution. In the linear velocity field approximation (34[ , the velocity field is a linear function of 
position, i.e. v(r) = Mr where M is a matrix. In general, M contains both anti-symmetric 
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TABLE II: N^/m 2 , in units of 10~ 4 GeW 2 for m x = 100 GeV 

Channel I II III IV 

WW, ZZ(0.73, 7.76) 38 30.8 7.9 0.6 

65(1,10.7) 51.9 41.8 10 0.5 

ft(l. 1,15.1) 57 45.4 9.8 0.3 

tm(0.95,6.5) 49.4 40.3 10.7 1.0 



TABLE III: JV 7 /m* in units of 10~ 4 GeV -2 for m x = 200 GeV 
Channel I II III IV 

WW, ZZ(0.73, 7.76) 13.5 11.0 3.1 0.4 
66(1,10.7) 18.4 15.0 4.10.4 

1,15.1) 20.2 16.4 4.3 0.3 

uu(0.95,6.5) 17.5 14.4 4.2 0.6 



and symmetric parts. If M is dominated by the anti-symmetric part (rotational flow), the 
infall is greatly simplified and the resulting caustics have the appearance of rings [3211 . I f M 
contains a significant symmetric part, the caustic geometry is more complicated. See [34J for 
a description. In this article, we assume that the caustics have the appearance of rings. If the 
flow has axial symmetry, the caustic ring is circular with constant cross section. Otherwise, 
the cross section will vary along the ring and the ring will not be circular. 

Let us consider cylindrical co-ordinates (s, z) where s = a/x 2 + y 2 . We assume axially 
symmetric infall about the z axis and reflection symmetry about the z = plane. We can 
then obtain an analytic solution for the dark matter density at points close to the caustic. 
With these assumptions, the caustic is circular with a tricusp cross-section described by the 
curve 



where 



(8) 



£ = 1 + 8^. (9) 
V 



Here, a is the caustic radius (s co-ordinate of the point of closest approach of the particles 
in the z = plane), p and q are the horizontal and vertical extents of the tricusp. See [32| 
for a detailed description. 

Let us use the self-similar infall model 51, 52| to estimate the mass infall rate for the 
flow of particles forming the caustic 

d 2 M v 2 
a ivi v t 

dndi = vf ^G (10) 

where v is the velocity of the particles forming the flow, / is a dimensionless quantity that 
characterizes the density of the flow and v rot is the rotation speed of the galaxy. 
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Let us define the two dimensionless co-ordinates R — (s — a)/p and Z = z/q. In the 
Z = plane, the density is given by 



p{R,0) ~ 0.34 
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Sfcp C and pfcp C are distances measured in kpc. For R ■ 
i? -1 / 2 for R > 0. For i? — > 1, the divergence ~ \(R 



the form 53 



• the dark matter density diverges as 
When Z ^ 0, the density takes 



p{R,Z) « 0.17 
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where Tj are the real roots of the quartic 
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The above formulae are only valid at points close to the caustic. 

The emission measure is calculated by integrating the density squared along the line of 
sight. Let b = n/2 — 9 be the galactic latitude. I is the galactic longitude chosen so that the 
galactic center is located in the direction I = 0, b = 0. We will assume that the caustics are 
spread over a distance ~ 10~ 4 pc. / is set equal to 2 x 1CT 2 [15]. The cut-off density close 
to the fold surface(near R = 0, Z = 0) is then « 2.15 x 10 3 /a v /p Gev/cm 3 . The density 
close to the cusp will be larger than this (the density falls off as the inverse of the distance 
to the cusp), but we will use the density close to the fold surface to set the density cut-off. 

Fig. [3] shows the emission measure averaged over a solid angle 10~ 5 sr for three different 
sets of caustic parameters, as a function of longitude I. b is set equal to and we assume that 
the caustic lies in the galactic plane. Figures 3(a), 3(b) and 3(c) are plotted for (a,p,q) = 



(7.5, 0.5, 0.5), (8.0, 0.1, 0.2) and (8.0, 0.1, 0.5) respectively with all distances measured in kpc. 
The earth's location is set equal to 8.5 kpc from the center. We expect the signal to be 
strongest when the line of sight is tangent to the ring. From the figures, we see that the 
emission measure is sensitive to the caustic geometry. The prominent features are the pair 
of peaks, or 'hot spots' separated by a few degrees. The first peak occurs when the line of 
sight is tangent to the fold surface (when s = a). The second peak occurs when the line of 
sight is tangent to the cusp line (when s = a + p). (In the limit p, q — > 0, the two peaks 
coincide, see [l2|). For the case when p = 0.5 kpc (Fig. 3(a)), the cut-off density was set 



equal to 400 GeV/cm 3 everywhere, while for the plots with p = 0.1 kpc (Figs 3(b) and 3(c) ), 
the cut-off density was set equal to 800 GeV/cm 3 everywhere. The magnitudes of < EM > 
for the hot spots depend on the values of the caustic parameters and also on the averaging 
scale (here chosen to be 10 -5 sr). Table IV shows the annihilation flux for the two hot spots 
for different values of the averaging scale Afl, for the case (a = 8.0, p = 0.1, q = 0.2) kpc. 
It is worth pointing out that if the triangular feauture in the IRAS map is interpreted as 
the imprint of the nearest caustic on the surrounding gas as in 44J, the implied caustic 



parameters are close to what we have assumed for Fig 3(b) 
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FIG. 3: Figures 3(a)[|3(b) 3(c) show the emission measure averaged over a solid angle A£l = 10 
sr for the three sets of caustic parameters. 



TABLE IV: Peaks of < EM > for (a = 8.0, p = 0.1,q = 0.2) kpc 
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(sr) 


(GeV/cc) 2 kpc 


(GeV/cc) 2 kpc 


io- 7 


1549.8 


2177.6 


10~ 6 


1048.3 


1283.3 


10" 5 


469.9 


718.1 


io- 4 


269.6 


239.4 


10" 3 


115.4 


52.0 



III. COMPARING THE SIGNAL WITH THE BACKGROUND 



The annihilation flux from caustics is thus given by 

< EM > 



S x 
110- 
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(1 



Let us compare this flux with the expected background. The EGRET measured background 



flux $>b g (Ei, E-z) from energies E\ to E 2 is given by [49j, |50| 

® bg {E 1: E 2 , 9, 0) = N y>bg (E u E 2 ) x N (9, 0) x 1(T 6 cm" 2 s" 1 sr" 1 (15) 

where 

^w^)=/j>(^)" 7 (i6) 



The function N (9,(f>) is energy independent and follows the fitting form given in |49J. For 
the four energy bands we have considered, N^^g = 198.8 for Band I (30 MeV - 100 MeV), 
28.9 for Band II (100 MeV - 1 GeV), 0.58 for Band III (1 GeV - 10 GeV) and 0.012 for Band 
IV (above 10 GeV). 

We now compare the caustic signal with the expected background. Since the gamma 
ray background falls off with energy (E~ 2 - 7 ) faster than the annihilation signal (E^ 1 - 5 ), we 
expect that the best chance for detection is at moderately high energies. At low energies, 
the background flux overwhelms the signal, while at very high energies, the signal is weak. 
We choose m x = 50 GeV since this choice gives the largest flux. For the quantity iV 7 /m 2 , 
we use the average value for the band. The averaging scale AQ is set to 10 -5 sr. 



Figures 4(a) and 4(b) show the expected annihilation flux (number of photons per square 
meter, per steradian, per year) as a function of angle I near the plane of the galaxy (b = 0) for 
the three sets of caustic parameters we considered, for Energy Bands III and IV respectively. 
This is contrasted with the expected diffuse gamma ray background. The figure shows only 
the flux from the caustic and does not include the smooth component of the halo. For a 
standard isothermal with core type of halo profile, we expect the signal from the smooth 
component of the halo to be much smaller than the signal from the caustic. However, if 
there are dark matter clumps along the line of sight, annihilation from the clumps may be 
significant In principle, the peaks in the signal and the sharp fall-off of flux are helpful 
in identifying the annihilation signal, particularly for the more optimistic caustic parameters 
and for small WIMP masses. For large WIMP masses, the annihilation signal is significantly 
smaller. 

Let us estimate the significance of detection by assuming an angular resolution AQ = 
10~ 5 sr, an integration time of 1 year, a detector area of 1 m 2 and an optimistic value 
of m x = 50 GeV. Further let us consider the 1 GeV - 10 GeV band, for which the mean 
value of ( m /iooGey) 2 * s ^1-1 (Tabled]). The number of photons received, assuming 100% 
detector efficiency, from a region with angular extent AQ = 10~ 5 sr is then (Eq. [T4l) . 
N 1>S ignai ~ 0.02 ( G ^j^2 kpc ■ Since an angular size of 10~ 5 sr corresponds to a spot of diameter 
roughly 0.2° at the equator, we may expect about 10 distinguishable spots in the angular 
range -0.1° < b < 0.1°, 68.2° < I < 70°. Here b is the latitude (vr/2 - polar angle 9) and 
I is the longitude. For the moderately optimistic case of caustic parameters (a = 8,p = 
0.1, q = 0.2) kpc, we find approximately 60 signal events and for the very optimistic case 
(a = 8,p = 0.1, q = 0.5) kpc, we find approximately 101 signal events. In this range, we 
find 716 background events. Assuming that the fluctuation in the background varies as 
the square root of the number of background events (N 7) f, g ), we may define the detection 

significance a = N JtS i gna i/ JNy~bg- For the case (a — 8, p — 0.1, q = 0.2) kpc, we find a = 2.2 
while for (a = 8,p = 0.1, q = 0.5) kpc, we find a = 3.8. A similar calculation shows 
that in the range —0.1° < b < 0.1°, 70.2° < / < 72°, we find 93 signal events for the case 
(a = 8,p = 0.1, q = 0.2) kpc, and 148 signal events for the case (a = 8,p = 0.1, q = 0.5) 
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kpc, with 700 background events. The values of a for the two sets of caustic parameters is 
found to be 3.5 and 5.6 respectively. 




I (degrees) / (degrees) 

(a) Energy Band III (b) Energy Band IV 

FIG. 4: (color online) The annihilation flux for the three different sets of caustic parameters 
(m x = 50 GeV), compared with the EGRET measured diffuse background. The background 
signal appears as a line of nearly constant slope. The characteristic peaks and the sharp fall-off 
of the signal may help in signal detection. Among the three curves that show the signal flux, 
the largest flux curve corresponds to the case (a = 8.0, p = 0.1, q = 0.5), the middle curve is for 
(a = 8.0, p = 0.1, q = 0.2) and the smallest flux is for (a = 7.5, p = 0.5, q = 0.5). Energy Band III 
is favored to Energy Band IV owing to the very small flux in Band IV. 



IV. CONCLUSIONS 

We calculated the gamma ray annihilation signal from a nearby dark matter caustic 
having the geometry of a ring with a tricusp cross section near the plane of the galaxy, in 
different energy bands. For such a caustic, the annihilation signal has two peaks, separated 
by a few degrees, depending on the size of the caustic. There is an abrupt fall-off of flux after 
the second peak. Since the diffuse gamma ray background flux falls off with energy faster 
than the signal, it is advantageous to look for the signal at moderately high energies. We 
compared the expected annihilation flux with the expected diffuse gamma ray background. 
The characteristics of the annihilation flux can in principle, be used to discriminate between 
the signal and the background. In practice however, we expect this to be a challenging task. 
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